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Exotica ex nihilo: Traversable wormholes & non-singular black holes 
from the vacuum of quadratic gravity 
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We present new traversable wormhole and non-singular black hole solutions in pure, scale-free B? 
gravity. These exotic solutions require no null energy condition violating or “exotic” matter and 
are supported only by the vacuum of the theory. It is well known that f{R) theories of gravity 
may be recast as dual theories in the Einstein frame. The solutions we present are found when the 
conformal transformation required to move to the dual frame is singular. For quadratic gravity, 
the required conformal factor is identically zero for spacetimes with = 0. Solutions in this case 
are argued to arise in the strong coupling limit of General Relativity. 


I. INTRODUCTION 


Recently, an exploration of the vacuum solutions of 
pure R^ gravity uncovered new black hole solutions, re¬ 
sulting, in part, from the lack of a generalized Birkhoff 
theorem [1]. As a theory, quadratic gravity has sev¬ 
eral intriguing features: it possesses scale invariance and 
with the addition of the Weyl curvature tensor = 
, is renormalizable [2]. Numerically it has 
been shown that the most general quadratic gravitational 
action, — aC^ for some constant a, also admits new 
black hole solutions [3]. Pure gravity is the only pure 
quadratic gravity that is ghost free HE] and, as an exam¬ 
ple of an f{R) theory, is sometimes dual to the Einstein- 
Hilbert action of General Relativity (GR) minimally cou¬ 
pled to a scalar field. We briefly review how this duality 
arises in section nn and discuss the cases where the dual¬ 
ity fails. It is precisely in this regime, where R = 0, that 
the new black hole solutions of [T] were found. The au¬ 
thors interpreted these solutions as a part of the strong 
coupling limit of GR. In this paper, we show that the vac¬ 
uum space of R^ gravity permits new traversable worm- 
hole and non-singular black hole solutions. We emphasize 
that these structures are supported only by the vacuum 
and do not require any unusual states of matter. 

In Einstein’s theory of General Relativity, every metric 
is a solution to Einstein’s equations for some associ¬ 
ated stress energy tensor T^y . It is therefore a challenging 
task to determine which solutions should be considered 
physical solutions. This dilemma lead to the develop¬ 
ment of ad hoc energy conditions intended to reasonably 
restrict properties of the sourcing matter. The weakest 
such condition, the Null Energy Condition (NEC), stipu¬ 
lates that the stress energy tensor of matter should satisfy 
T^yk^k^ > 0 for any arbitrary null vector k^. Impos¬ 
ing such conditions prohibits the construction of many 
spacetimes in the context of GR, including traversable 
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wormholes [6]. Wormholes have appeared in many sci¬ 
ence fiction settings, recently giving rise to additional 
research on the subject [7]. Besides a means of rapid 
interstellar travel used by advanced civilizations, worm- 
holes are discussed in the quest to understand the rela¬ 
tion between entanglement and the possible emergence of 
spacetime [8] . Such ideas were applied to the information 
paradox, yielding the so called, ER = ERR conjecture 
[9] , which suggests that entangled particles are connected 
via a (non-traversable) Einstein-Rosen bridge. Likewise, 
it has long been speculated that quantum gravity may 
somehow resolve singularities in black hole spacetimes 
leading some researchers to the idea that the Universe 
itself, might have been created on the interior of a black 
hole (for some early work see pI)U2Q] L 

In section [HI we discuss how the new spacetime solu¬ 
tions are possible and present our solutions in section [ITU 
We then show explicitly in section IV that observers and 
light rays can traverse the throat of the wormhole and 
discuss some of the trajectory’s properties in section |V| 
Einally, we discuss the appearance of the wormhole to 


asymptotic observers in section VI 


II. GRAVITY AND ITS VACUUM 
STRUCTURE 

The R^ gravity action is a specific example of the more 
general f{R) action [27], 

r 

S = ^ J + Sm- (1) 

When the matter action Sm exhibits conformal symme¬ 
try and E{R) = df /dR is non-vanishing, it is well known 
that this action is equivalent to Einstein gravity mini¬ 
mally coupled to a scalar field. This is shown by writing 
the action in terms of the metric g^y = D?g^y, where 
the conformal factor is given by U = F(R)^/^. There¬ 
fore, any solutions found in such f{R) theories will also 
be found in Einstein gravity coupled to a scalar field. 
However, as was shown in [1], this is not necessarily the 
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case when F{R) vanishes, and the metric is ill de¬ 
fined. Wormhole solutions in f{R) have been studied 
extensively [2T1426] . 

In f{R) gravity the metric has equation of motion, 

F{R)R^, - - V^V,F(i?) 

+g^,OFiR) = M-%,. ( 2 ) 

In the above, is the stress-energy tensor of any 

external matter present Q Focusing on vacuum solu¬ 
tions, = 0, we notice that any metric producing 
f{R) = F{R) = 0 will be a solution of this theory. These 
conditions may constraint R but leave significant freedom 
in the form of the Ricci tensor R^y. 

Focusing on the simplest case which satisfies the above 
criteria, f{R) = we find that every spacetime with 
traceless Ricci tensor is a solution to the R^ gravity 
EOM. This abundant freedom in the form of R^y is how 
Birkhoff’s theorem is circumvented. In standard Einstein 
gravity, R = 0 implies R^y oc T^y and therefore any met¬ 
ric sourced by matter with traceless energy-momentum 
tensor will also be a solution to the vacuum of R^ grav¬ 
ity. One notable example is the electromagnetic field with 
Fm = —\F^yF^^. We will show that this new freedom 
found in the vacuum of R^ gravity allows bundles of ingo¬ 
ing modes to evolve from converging to diverging modes 
in section |V| In ordinary Einstein gravity, such evolution 
would only be possible through a violation of the NEC 
as can be shown through the Raychaudhuri equation-an 
equation which describes the evolution of the bundle’s 
divergence 0. Eor null geodesics with tangent null vector 

^ the Raychaudhuri equation is 

(3) 

Here a^y is the shear tensor and uj^y the vorticity ten¬ 
sor. The vorticity can always be set to zero by choos¬ 
ing a coordinate system in which the congruences are 
hypersurface orthogonal. The remaining terms on the 
RHS are all negative in GR if the NEC is satisfied 
{R^yk^k^ = Mp‘^T^yk^k^ > 0). This implies that 
< 0 and hence bundles of null geodesics are always 
converging. In R? gravity, vacuum solutions can have 
R^yk^k^ > 0 so that the above conclusion need not ap¬ 
ply- 

III. NON-SINGULAR BLACK HOLES AND A 
TRAVERSABLE WORMHOLE 

The most general static and spherically symmetric 
spacetime is given by 

ds‘^ = —fi{r)dt‘^ + f 2 {r)dr‘^ + (4) 


^ We use the reduced Planck mass, = I/SttGn, metric signa¬ 
ture ( —, +, +, +), and □ = V^V^. 


where dU^ = dO^ + sin^6>d0^. We make the following 
metric ansatz 

/i(r)=G(r) and, f2{r) = ^2 1 fc2 ( 5 ) 

Introducing a reparametrization of the radial coordinate, 
p _ ^2 _ ^ 2 ^ for some constant /c, the metric becomes. 


ds‘^ = -G{l)dt‘^ + Tva (6) 

G[l) 


If G{1) 1 as / ^ ± 00 , this metric connects two asymp¬ 

totically Minkowski spacetimes. The constant k will set 
the minimal wormhole throat radius as r cannot become 
smaller than this value. 

The Ricci scalar of this spacetime is given (in terms of 
G)by, 

- + l^fR = i((P + ffGp 

+{k^ + 1^){G - 1) + Gk^ , (7) 

where ' = d/dl. Requiring R = 0, yields a two parameter 
solution for G. A complicated exact solution of this sec¬ 
ond order ODE is given in Appendix [A| Eor our purposes 
we need only discuss certain limiting behaviors. Taking 
the large I limit of Eq. 0, corresponding to large r, and 
solving the R = 0 ODE, we find the asymptotic solution. 


^ ±(^) d - ^ + (8) 

This resembles a Reissner-Nordstrom metric of mass M± 
and charge Q±. These need not be the same for the two 
different asymptotic regions. Hence we label them by a 
± to denote the values measured by observers located 
at I Too. Note that these parameters arises as in¬ 
tegration constants and should not be associated with a 
conserved “mass” and “charge” even though we denote 
them by M and Q. Eor example, the approximate solu¬ 
tion in equation 0 drops a non-analytic term, seen only 
in the large I limit of the exact solution, that goes as 
^2iogdA) ^ Hence taken as the coefficient of the 1/P 
term in the series expansion, cannot be a well defined fi¬ 
nite charge. We will discuss this point further in section 

ED 

At small /, the solution is approximately. 


G{1 « 0) i [l + 2(Go — 0.5) cos{2l/k) + kvo sm{2l/k)] 

^ (9) 

where Go = G(0) and vq = G'(0) are the initial condi¬ 
tions determining the shape of the throat. The constants 
M± and Q± are expressed in terms of these parameters 
in section |VI[ Given appropriate choices of Go, vq^ we 
ensure that G is strictly positive for all 1. The solutions 
are shown in figure 0 . In these cases the spacetime has 
no horizons even though it appears as a black hole to 
asymptotic observers. 
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For certain initial conditions our solutions can develop 
horizons, however the solutions remain free of singular¬ 
ities. To the asymptotic observer these solutions ap¬ 
pear as non-singular black holes (or one-way traversable 
wormholes). In general there can be up to two regions 
where the I = const slices become spacelike, whenever 
G becomes negative. This can be seen by noting that 
for \l/k\ > 1 the ODE is similar to an overdamped oscil¬ 
lator with equilibrium at G = 1 and therefore any zero 
crossings must occur in the range I G [—k^k]. As 0(1) 
changes in G occur on a length scale / ^ /c, we expect 
G to oscillate at most once during this interval, possibly 
below zero, and then it will asymptote monotonically to 
its equilibrium value. This may occur on either side of 
the I = 0 point and there are at most two regions where 
G is negative. Penrose diagrams show the causal struc¬ 
tures of the three spacetimes discussed above in Figure 
[21 In the remainder of this work we focus on the cases 
with no horizons, where G remains well behaved so that 
the {t, / , ^, 0} coordinate system covers the entire space- 
time manifold. In other words we will assume that G is 
a continuous non vanishing function of I that asymptotes 
to 1 as / ^ Too. 



FIG. 1: Plot of G{l/k) vs l/k. The blue/dashed, 
red/dotted-dashed, and green/solid have initial 
conditions (Go,'Co) set as (1,0), (0.2,0),and (l,l//c) 
respectively. 


IV. GEODESICS OF TIMELIKE AND NULL 
OBSERVER 

Here we compute the paths of null and timelike 
geodesics traveling in the radial direction {dft = 0). We 
show that a timelike observer can, in principle, safely 
traverse the wormhole. We begin with a discussion of 
null geodesic which have ds‘^ = 0 and are therefore 
parametrized by, 

ft = 

Assuming no horizons, this may be integrated to find 





FIG. 2: Penrose diagrams of the three possible causal 
structure found in the vacuum solutions of ansatz (H- 
The shaded regions have G < 0 and the I = const slices 
are spacelike. Top: G is well behaved and no horizons 
exist. Middle: G is negative in a single region between 
I G [/_,/+]. Bottom: G is negative in two regions 
between I G and I G 
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the path l{t) of radially propagating photons. In sec¬ 
tion [V| we shall see that the affine parameter defined by 
dt/dX = cG~^ for some arbitrary constant c is a mono¬ 
tonic boundless function of t. Hence the light rays tra¬ 
jectory through the throat can be described in a finite 
range of the affine parameter. 

We now turn to timelike geodesics. Given a four ve¬ 
locity 14^ (r) with r being the observer’s proper time the 
geodesic equation reads, 

= 0 . ( 11 ) 


The normalization of u for radial trajectories gives the 
following condition, 


r dt\ 

i'- ^ 1 


\dT J 

1 Q2 1 

Kdr) 


1 

+ G- 


( 12 ) 


Eliminating = dt/dr from the I component of the 
geodesic equation, one can solve for /(r) through. 


dr^ 2 


(13) 


and t{r) is determined from Eq.(12). To visualize the 
geodesic motion, we embed the geometry in a higher di¬ 
mensional space. Suppressing the 0 coordinate and tak¬ 
ing a time slice at a particular t (which looks the same for 
any choice of t), we embed the wormhole in a 3D space 
as shown in figure (3a). This particular solution corre¬ 
sponds to a Go = 1 and t’o = 0 wormhole. The red line 
represents the radial trajectory of a timelike observer. 
We also show how different values of Go, vq (defined and 
color coded in the caption of figure 0 ) affect our coura¬ 
geous explorer along his journey. 

Note that the solutions can have regions where grav¬ 
ity is attractive or repulsive towards / = 0. The tran¬ 
sition occurs wherever G' switches sign. Eor instance, 
our solutions with Go = 1 have repulsive gravitational 
forces near the throat and incoming observers must pos¬ 
sess sufficient velocity to traverse the throat. In figure 
(3b), our infalling explorer labelled by the solid green 


curve does not have sufficient speed and is repelled out 
of the wormhole. However, the solution with Go = 0.2 is 
attractive everywhere towards the wormhole and the ob¬ 
server travelling in this background (designated by the 
red dot-dashed line in figure ( |3b| )) oscillates back and 
forth through the wormhole. It is even easy to find so¬ 
lutions where G asymptotically relaxes back to 1 from 
above. In these cases our asymptotic observers would 
associate a negative M± to the wormhole and feel a re¬ 
pulsive force from it! This might be concerning at first 
sight as this would imply a spacetime with negative mass 
in standard GR. However the situation is more compli¬ 
cated in gravity as we will discuss in section VI 


V. TIDAL FORCES THROUGH THE 
WORMHOLE 

In this section we determine how bundles of geodesics 
evolve as they pass through the wormhole. The non-zero 



(a) Embedding of the wormhole in higher 
dimensions with a throat-traversing, radial 
trajectory traced in red. We have 
suppressed the 0 coordinate such that every 
circle at some constant z is a sphere having 
that circle’s radius. 

//k 




(b) The coordinates I and t as seen by observers 
with identical initial speeds at / = 0 and 
attempting to cross the different throats. The 
observer on the green geodesic fails to traverse the 
wormhole and is spit back out. The observer 
labelled by the red dot-dashed line oscillates 
around the everywhere-attractive wormhole. The 
blue dashed observer successfully passes through 
the throat. 

EIG. 3: Wormhole structure. 
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components of the Ricci tensor read, 
21G' 

r;( - 

2 


Rt 


2^[wtp+^ j’ 

( 2 G)-i( 




Rgg = l-G- IG', 


k^ + P 

_^f4k‘^G + 2l{k‘^ + P)G' 

(fc 2 + ; 2)2 ^ 

= sm{9)‘^Rgg. 


G" 


(14) 

(15) 

(16) 


Taking the radial geodesics with affine parameter A, 
k = ^( 1 ,±G, 0 , 0 ) for in(-) and outgoing(+) geodesics 
respectively, we find 




/af \2 2k^G^ 

V^J {iPTW' 


(17) 


This is strictly negative and therefore can allow for 
> 0. We now compute the divergence explicitly. The 
divergence of a bundle of geodesics with tangent vectors 
is given hy 0 = Vak^ = . To deter¬ 

mine k^ = §y one must look at the geodesic equation for 
the null ray. The k^ equation is, 




— fc* + — fc* — 0 

dr + G^ 


(18) 

(19) 


where in the second line we assumed ^k^ = 0 , appro¬ 
priate for a time independent metric. This is solved by 
k^ = cG~^ for some integration constant c which denotes 
the arbitrary choice in the affine parameter’s normaliza¬ 
tion. We take the value c = 1 . In the cases were G 
is well behaved so that our coordinate system is valid 
everywhere, we can compute 0 for any /, 

The ± corresponds to in and outgoing modes respec¬ 
tively. As expected, 0 switches sign when the light rays 
cross the throat at / = 0. Moreover, since 0 gives the ex¬ 
pansion/contraction of the cross sectional area of a bun¬ 
dle of light rays (it is an area, as orthogonal hypersur¬ 
faces to null geodesics are two dimensional), the radial 
null geodesics of the metric in Eq. ^ are not affected 
by the form of G. 

A timelike observer will also feel increasingly squeezed 
by the tidal forces as we enters the wormhole. Unlike 
for null geodesics, such observers will be sensitive to G. 
Morris & Thorne [ 6 ] used these tidal forces to impose 
constraints on the wormhole spacetime by requiring that 
a human sized observer feels less than a ^ of force. We 
avoid imposing any such conditions here but give a brief 
derivation of the tidal forces below. We refer the reader 
to [ 6 ] for additional details. 

Consider an observer O' at rest with respect to 
{/, 6 >, 0 }, and proper reference frame with basis vectors. 




1 




ei = 


sm 


Oy/k^ + 


( 21 ) 


In this coordinate system the metric takes the simple 
form = (— 1 , 1 , 1 , 1 ). A second observer O" located 
at the same point, but moving with a speed v in the 
radial direction, will also have his own reference frame 
whose basis can be written in terms of the coordinates of 
O' through a Lorentz boost, leaving 77 ^^ unchanged. 


^6 = + 7 ^ 67 , 

ei = yve^ + ye^. 

Here the gamma factor is the usual 7 = (1 — 7 ;^)“^/^. 
In this frame the moving observer has four velocity u = 
Gq. Since the four acceleration is perpendicular to the 
four velocity, our radially moving observer’s acceleration 
will only have a non zero radial component in his own 
reference frame, namely a = aei . To find a is to write a = 
in the {t,/, 6 >,(/)} coordinates which have constant 
vielbeins and use. 


oiUt = a • Gt = aGi • Gt. 
Here, Ut = and. 


a = — 


1 



d'f' 


(23) 


(24) 


At I ±00 and assuming 7 1 , the acceleration felt 

is a ~ For I « 0, the term G' contains terms 

evolving as G^/k and vq. Hence one may tune these 
values so that the observer crossing the throat will feel a 
desired acceleration. 

The tidal acceleration felt by the observer O" can also 
be computed in his reference frame. The relative accel¬ 
eration of two nearby points separated by ^ = ( 0 , ^) as 
measured by O" is given by. 




(25) 


with, 


G" 

-^ioio ~ ^ ’ -^2626 “ -^3036’ 


-^3036 “ R 


IG' 


3030 ’ P + k^ 


2 2 / 1(1 +r)G' + 2PG 


■ 7 V 


V 2(^2 + fc 2)2 


)• (26) 


We can explicitly see how the size of the throat k influ¬ 
ences the strength at which our traveler is squeezed in 
the angular directions. 


VI. THE WORMHOLE TO ASYMPTOTIC 
OBSERVERS 

From Figure Q we see that there is a lack of parity 
symmetry in the I component when 7 ;o 7 ^ 0 . As previously 
mentioned, this causes observers in different asymptotic 
regions to see different looking wormholes, characterized 
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by M± and Q± seen at I = zboo. These are determined 
by Go, '^0 from expanding the exact solution of G and 
picking out the asymptotic coefficients. Defining u = 1//, 
we expand the solution of G around = 0 and define, 


M± = lim (27) 

u^±o z au 

d2 

Q±= lim ^G(u,Go,vo). (28) 

u^±o au^ 


One finds that, 


=F2M± Mo sin(V37r/2) 

—-— = “2— + cos{V37r/2)Go ±- ^ - kvq 

(29) 

with Mq ~ 0.7 being a number whose exact expression 
is given in Appendix 0. Note that if G is to approach 
one from below in the I ^ —oo limit, then the sign of 
M_ will be negative but this is just because our ’’radial” 
coordinate I takes on negative values. The ± sign in front 
of vq is responsible for the asymmetry. Even though this 
implies that asymptotic observers will feel different ac¬ 
celerations toward the wormhole, one cannot gain energy 
by going from one region to the other. The asymmetrical 
shape of the throat (as is seen by the solid green curve 
in FigQ will prohibit this as can be shown by Eq.(13) 
which implies (dZ/dr)^ + G is a conserved quantity and 
so Av = [dl/dT]i=±oo = 0. 

Meanwhile the coefficient is found to be. 


= -^±v^sin(A/37r/2)Go—cos(\/37r/2)/ci;o—21og('uk). 

(30) 

The very last term has a logarithmic dependence so 
that, not only is not a constant, it diverges as 
14 ^ 0. This shows that the exact solution and the Reiss- 
ner-Nordstrom spacetime differs radically in the asymp¬ 
totic region. In the case of GR with electromagnetism, 
the conserved charge is defined by the conserved current 
integrated over a spacelike hypersurface. 
Here no such analogue electromagnetic field is present. 
The logarithm found is just a feature of our particular 
wormhole solution. Choosing a different ansatz, such as 
simply replacing the angular part of Eq. by, 

= (l^ + (31) 

allows us to obtain wormhole solutions that asymptoti¬ 
cally approach Reissner-Nordstrom metrics without log¬ 
arithmic contribution. 

Conserved charges are tightly linked with symmetries 
of the theory. In our case, we have a timelike killing 
vector Bt = -^ with which we can construct a conserved 
current and find the corresponding charge; however, the 


interpretation of charges in gravity is tricky. Pre¬ 
vious attempts [28l |29] by Deser & Tekin, to define an 
energy for higher curvature gravity in analogy with GR 
have shown that every asymptotically flat solution of 
gravity has vanishing energy. An alternative definition 
was proposed by the same authors and relied instead on 
the corresponding Poisson equation of higher curvature 
gravity m- In terms of this second definition, the lead¬ 
ing term of the source in the Poisson equation defines 
the energy and allows for a less degenerate classification 
of spacetimes but still yields a vanishing energy for solu¬ 
tions with everywhere vanishing Ricci scalar. Therefore 
even though an asymptotic observer will feel that the 
wormhole is exerting a gravitational force (just as if it 
were a spherical object of mass M± in ordinary GR) he 
shouldn’t associate this value with a conserved charge 
such as the energy. Hence one might wonder if these 
quantities can evolve upon the introduction of small per¬ 
turbations. 


VII. CONCLUSION 

In this paper we have presented new static wormhole 
and (non-singular) black hole solutions found in the vac¬ 
uum of R^ gravity. The traversable wormholes are are 
solutions in D = 4 which do not require the support 
of any NEC violating exotic matter. Their existence 
could prove important as it is argued in [1] that the new 
R = 0 solutions might be part of the strong coupling limit 
{Mp 0) of General Relativity. It would be interesting 
to study the role of such spacetimes in this limit using 
other tools. For instance it was shown recently that in 
the context of AdS/CFT, the strong subadditivity condi¬ 
tion of entanglement entropy imposes constraints on the 
bulk geometry of certain spacetimes [31]. One would re¬ 
quire a dual GET to asymptotically R = 0 solutions to 
apply similar results. 

In the context of our particular solution, we have 
shown that asymptotic observers view the wormhole as 
a black hole with “mass” M± and varying “charge” Q±. 
We have stressed that one should not be associating these 
values as conserved charge. They are taken as coefficients 
of the 1/r and 1/r^ terms in the metric expansion and 
do not correspond to some conserved quantity as was the 
case in GR. This raises the question of whether M± and 
Qi, along with the whole wormhole itself, are dynamical 
once small perturbations are introduced. It would also 
be interesting to look for further exotic solutions to gen¬ 
eral f{R) theories which cannot be captured by Einstein 
gravity coupled to a scalar field. 
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Appendix A 


The general solution to G is given by 


Q(^l^ ^ ^/c COS arctan(///c)) ^ ^/c sin (v^arctan(///c)) 


+ P 


V3Vk^ + P 


kvo 


6Vk‘^ + - 


- 2 V 3 - 3 + VSTTtan ((1 + 

+ (_3+v^) - V3^fi(-i+V3) + 


2k 

{3 + V3){k-il) 


11 ‘Ik 

(2 + a/s) 2 -f’i(l, 2(1 “ 2^^ ~ ^ ~ k-il^ 


1 1 

+ 2-F’i(1, 2(1 + ^ ~ k -iV 


(Al) 


Here 2-^1 («, b; c; is the ordinary hypergeometric function. The coefficients M± and Q± can be read from the series 
expansion of G at large I as written in Eq. (27, 28). The result is, 

T2M± ^ sin(\/37r/2) 

= cos(v37r/2)Go ±-7=-fcuo 


fc 


73 


+ 


273 


cos 


,(737r/2) + 273 + 4 - TTtan ((1 + 73)7r/4) , {A2) 


^ = ± 73sin(737r/2)Go — cos(737r/2)fcuo 


fc2 


+ - sin(737r/2) + 273 + 4 - Trtan ((1 + 73)7r/4) 

- 21og(2uk) - 2(1 - je-m) - ((1 - 73)/2) - ((1 + 73)/2) 


(A3) 


Here ^e-m is the Euler-Mascheroni constant. The functions and are Harmonic numbers and polygamma 

functions respectively, both are defined from the Gamma function through, 


Hn = 'JE- 


M 


r'(n + l) 
r(n + l) ’ 

= r'(7/r(7- 


(A4) 

(A5) 
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